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Abstract 

The thermal entanglement of the Hubbard dimer (two-site Hubbard model) has been studied 
with the nonextensive statistics. We have calculated the auto-correlation (Og), pair correlation 
(Lq), concurrence (Tq) and conditional entropy {Rq) as functions of entropic index q and the 
temperature T. The thermal entanglement is shown to considerably depend on the entropic index. 
For q < 1.0, the threshold temperature where Tq vanishes or Rq changes its sign is more increased 
and the entanglement may survive at higher temperatures than for q = 1.0. Relations among 
Lq, Tq and Rq are investigated. The physical meaning of the entropic index q is discussed with 
the microcanonical and superstatistical approaches. The nonextensive statistics is applied also to 
Heisenberg dimers. 
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I. INTRODUCTION 



Quantum entanglement is one of the most intriguing subjects in quantum information 
theory, and it has been investigated from various viewpoints in the last decade (for a review, 
see Refs. [l, 2|, related references therein). Quantum entanglement is expected to play an 
essential role as a resource in quantum communication and information processin g. M any 
studies have been made on quantum entanglement with quantum spin models jsl-jlll and 



fermion systems 12|- 27| to clarify both its qualitative and quantitative aspects. The in- 
terface between the quantum information and statistical mechanics has been considerably 
investigated in recent years. It has been shown that entanglement of two neighboring sites 
shows a sharp peak either near or at the critical point where the phase transition takes place 
sj-Q, Q . This sug gests an intimate relation between quantum entanglement and quantum 



phase transition 



nis sug gest s ai 
ition P, I9I-2I 



26|. 



These studies mentioned above have been made within the Boltzmann-Gibbs statistics 
(BGS). In the last several years, there is an increased interest in the nonextensive statis- 
tics (NES), which was initially proposed by TsaUis |28l-l30| . This is because the standard 
method based on the BGS cannot properly deal with nonextensive systems where the physi- 
cal quantities such as the energy and/or entropy of the iV-body system are not proportional 
to A^. The nonextensivity has been realized in the three typical systems: (1) systems with 
long-range interactions, (2) small-scale systems with fluctuations of temperatures or energy 
dissipations, and (3) multi-fractal systems [3l|- For example, in a gravitating system with 
;he long-range interaction, which is a typical case (1), the specific heat becomes negative 



32l |. A cluster of 147 sodium atoms, which belongs to the case (2), has been reported to 



show the negative specific heat 
entropy) proposed by TsaUis [28 1 



The generalized entropy (referred to as the TsaUis 
30| is given by, with Tr p = 1: 
Tr pi - 1 



Sa — k 



1-q 



(1) 



where A; is a positive constant, p the density matrix, Tr the trace and q the entropic index. 
In the limit of g — )■ 1, Eq. ([T]) reduces to the von Neumann form. 



5*1 = —k Tr plnp. 



(2) 



where Tr p = 1. We will set k = ks {ks'- the Boltzmann constant) when we discuss the 
thermodynamical entropy. The nonextensivity (non-additivity) in the TsaUis entropy is 



explained as follows. For a system consisting of two independent subsystems A and B with 
the density matrices, p{A) and p{B), the density matrix of the total system is described by 
p{A, B) = p{A) ® p(i?) from which the Tsallis entropy is given by 



S,{A,B) = S,{A) + S,{B) + 



k 



B 



SM)S,{B), 



(3) 



Sg{ri) denoting the entropy of the subsystem 77 (= A, B). Eq. ([3]) shows that the entropy is 
additive for g = 1 and non-additive for 1. The quantity (g — 1) expresses the measure of 
the nonextensivity. The NES has been applied to a wide class of subjects such as physics, 
chemistry, and biology 



34, 



35|. 



Small-scale systems belong to nonextensive systems as mentioned above. It is necessary to 
take into account the effect of the nonextensivity, when we study the properties of quantum 
information of qubits, which have been mainly investigated within the BGS [l2|-j27|. In 



recent years, the NES has been applied to a study on quantum information [36|-[53|]. It has 



been shown that the non-additive Tsallis entropy yields a better measure for the separabi 



36 



3 



40 



53|. Ref. 



ity 



39| 



criterion of entanglement than the additive von Neumann entropy 
discusses an NES generalization of the von Neumann mutual information which is shown 
to strongly depend on the entropic index q. The above list of applications of the NES to 
quantum information is not complete: relevant references are 



presented in 



34 



35 



It is worthwhile to apply the NES to the Hubbard model [5^ , which is a typical model for 
strongly correlated fermion systems and which has been widely adopted for a study on quan- 
tum information [l2|-j23]. Despite the simplicity of the model, however, it is very difficult 
to obtain its exact solution except for some limited cases. In order to obtain a reasonable 
solution, various approximate methods have been proposed and developed. The Hubbard 
model provides us with good qualitative description for many interesting phenomena such 
as electron correlation, magnetism, superconductivity and quantum entanglement. We may 
employ a finite-size Hubbard model to obtain an analytical, exact solution. The Hubbard 
dimer (two-site Hubbard model) has been adopted as a simple model which can be analyt- 
ically solved. Although the Hubbard dimer seems a toy model, it h as play ed an important 



iQei, It naspiay 
tion [14, His 



role as a model of qubits in the theory of quantum information [iJ, ll5|, Il8 |. 

In our previous papers SS-SsJ, we applied the NES to Hubbard dimers to investigate 
effects of the nonextensivity on their thermodynamical and magnetic properties, bearing 
small-size systems in mind. It is interesting to examine the effect of the nonextensivity on 



3 



the properties of the quantum entanglement of two qubits [2] described by the Hubbard dimer 
within the NES, which is the main purpose of the present paper. Among various quantities 
expressing thermal entanglement, we have calculated the pair correlation, concurrence 



591, 



GOj and conditional entropy 6l|, |62|]. As will be shown in our study, the entropic index q 



has considerable effects on the properties of thermal entanglement which may be improved 
by the nonextensivity. The concurrence of the Hubbard dimer has been discussed within 



the BGS 



14 



15|, 



proposed in Refs. 



18| ■ Th e generalization of the conditional entropy to the NES has been 
40, 



45( 1 ■ This is the first systematic study on the thermal entanglement 



of nonextensive fermion systems as far as we are aware of. 

The paper is organized as follows. In Sec. 2, after briefly reviewing the maximum-entropy 
method (MEM) in the NES 28|], we derive the density matrix to obtain auto- and pair 
correlations, concurrence and conditional entropy. In Sec. 3 we present model calculations 
of relevant quantities as functions of the entropic index q and the temperature T. In Sec. 
4, we make a comparison among the pair correlation, concurrence and conditional entropy. 
Effects of magnetic field and interatomic interactions in the adopted model are investigated. 
The physical mea ning of the entropic index q is discussed with the use of the microcanonical 
approach (MCA) |63l-l67| and superstatistical approach (SSA) [68|-[70|. Sec. 5 is devoted to 
our conclusion. 



II. FORMULATIONS 
A. Hubbard dimers 

We consider the extended Hubbard dimer whose Hamiltonian is given by 

2 

H = -t^{al^a2a + 4^ai^) + ^jt'^j; + XI '^^'^'^2,7 + ^ nia^2-<7 

2 

- /^B^XlKt -^iJ, (4) 

where rija = a^j^aja, dja denotes the annihilation operator of an electron with spin a (=t, |) 
on a site j (=1, 2), t the hopping integral, U the intraatomic interaction between electrons 
with opposite spin, Vi (V2) the interatomic Coulomb interaction between the same (opposite) 
spin, the Bohr magneton and B an applied magnetic field. By using the standard basis 

4 



for half-filled case with two electrons given by 

l^l) = I ti)l|0)2, 1^2) = |0)l| ti)2, 1^3) = I t)l| 1)2, 

\^,) = |i)l|t)2, \^,) = \t)l\t)2, \^e) = \i)l\i)2, 

we obtain the Hamiltonian matrix, 



H 





u 





-t 


-t 





\ 







u 


-t 


-t 










-t 


-t 


V2 













-t 


-t 





V2 



















Vi 







v 

















Vi + 2fiBB J 



Six eigenvalues of the system are given by [15|, l71| 



e, = liU + V2-D), liU + V2 + D), U, V2, - 2/is5, ^1 + 2^^5 

for i = 1 to 6, 



and the corresponding eigenvectors are given by 



$1) 
$2) 

$4) 

$6) 



sin^ 



(|ti)l|0)2+|0)l|ti)2) 



COS 6 

71" 



(|t)l|i)2+|i)l|t). 



V2 

(I ti)i|o)2 + |o)ii ti)2) - (I t)ii ;)2 + 1 1)2) 



V2 

^(in)i|o)2-|o)iin)2) 



V2 



(|t)lU)2-|i)l|t)2), 

lit). 

Ii)l|i)2, 



V2 

lt)l|t)2, 



where 



tan^ 



At 



D = ^{U - V2f + 16^2. 



For t/U < 1 with Vi = ^2 = 5 = 0, we obtain 
ei = 



sin^^ 



4^2 4^2 
— ^, 62 = t/ + — , 63 = f/, 64 = 65 = eg = 0, 

W o . 4t2 



t/2 



cos'0 



(5) 

(6) 

(7) 

(8) 

(9) 
(10) 

(11) 

(12) 
(13) 

(14) 
(15) 
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where ei is the lowest eigenstate for f/ > 0. 

The partition function in the BGS is given by 



^BG 



Zi = Tre-'^^ = ^ 



^ 2 e-/^(f^+^2)/2 cosh (^^^ + e-^^ + e"^^^ + 2 e"'^^^ cosh(2^/iB5), (16) 

where /3 = 1/kBT. From Eq. (JT6ll we can obtain various thermodynamical quantities of the 
system. 

B. Maximum-entropy method in the NES 



We will study the Hubbard dimer given by Eq. 



ID within the NES, where the PDF 
le Tsallis entropy. At the moment 



281], (b) un- normalized method |29| . 



or the density matrix is evaluated by the MEM for t 
there are four MEMs in the NES: (a) original method 

(c) normalized method jsoj], and (d) optimal Lagrange multiplier (OLM) method [7^. The 
four MEMs are compared in Ref. |3]|. Among the four MEMs, (c) normalized MEM and 

(d) OLM-MEM with the g-average have been widely adopted for a study of nonextensive 



systems, because they are be 



ieved to be more superior than (a) original MEM 28|] and 



(b) un-normalized MEM 29|, |31|- Recent papers 73N75| . however, have pointed out that 
thermodynamical quantities obtained by the g-average are unstable for a small change in 
the PDF, whereas those obtained by the normal average are stable 3]. The stability of the 
g-average is currently controversial [3]- [3]. Althou gh ( c) normalized MEM [s^ with the 



g-average was adopted in our previous papers [55 



58| |79[| . we have employed in the present 
study, (a) original MEM with the normal average 2^, 2^. In Appendix A, thermodynamical 
quantities of the entropy, specific heat and susceptibility calculated by (a) original MEM 



28 , |73| are summarized and compared to previous ones obtained by (c) normalized MEM 
30l | with the g-average 55l-l58l| 791]. In Appendix B the NES with (a) orig mal MEM 0, Q 



is applied also to Heisenberg dimers. 
Imposing the two constraints given by 



Tr (p, H) = {H), = E, 



(17) 
(18) 
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we obtain the density matrix given by 



with 



p, = — ExpJ-/3(/7-E,)], 



X, = TTExp\-(3iH-E,)], 



(19) 



(20) 



where {■)q expresses the average over pg, [3 the inverse of the temperature and Exp^(x) is 
defined by [73| 



Exp^(x) = [l + {l-l/q)x]f'-'\ (21) 

with = max(?/,0). Note that Expg(x) is diff'erent from the conventional g-exponential 
function expg(a;) defined by 



exp (x) = [1 + (1 - g) 



X 



1/(1-9) 



The two g-exponential functions, Expg(x) and expg(x), have the relation 73| : 
expg(a;) = Exp2_q((2 - g)a;), Expg(a;) = exp2_g(a;/g). 



(22) 



(23) 



Both Expg(x) and expg(x) reduce to the exponential function exp(x) in the limit of g — )■ 1.0. 



C. Auto- and pair correlations 

For the Hubbard dimer under consideration, we obtain 



Wi 



where 



Wi 



ExpJ-/3(e,-E,)]. 



(24) 

(25) 
(26) 

(27) 



The energy Eq in Eq. ( I25l) includes the partition function Xq which is expressed in terms 
of Eq in Eq. (!27|) . Then Eq and Xq are self-consistently determined by Eqs. (!25|) -( l27|) for 
given g and (3. 



We first consider auto- (Oq) and pair correlations (Lg) defined by 

2 

Lg = ^ {ni„n2a - ni„n2-a)g ■ 
When we employ the relations given by 

EKt'^^-^). = (If ) ' 

j=i \ / q 



Eqs. 1^^ and become 



dVi/ 



^ _/8H\ 



We may evaluate {dH/dO)^ with 6 = U , Vi, and 1^2 as follows. Taking the derivat 
Xq in Eq. ( l20l) with respect to 6*, we obtain 

<9Xg . „..,fdH dE, 



de 



/3Tr{(ExpJ-/3(i7-i^,)]) (f^-f^)}, 



which leads to 



/dH\ _ dEg^ 1 dXg 

x'W/^ ~ ~de~~ pXg de ■ 

From Eqs. fl25|) - fl271) . self-consistent equations for dEg/dO and dXgjdO are given by 

dEg _ dEg dXg 

dXg _ dEg dXg 

with 



8 



where an explicit expression for cie is not necessary (see below). Solving Eqs. ([38 
respect to dEqjdO and dXq/dO and substituting them to Eq. f l37|) . we obtain 



dH\ ^ _C2e_ ^ ±_ 

99 021 Xg 



El dei 



Simple calculations with the use of Eqs. (E]) and fHT]) lead to 



dH\ _ 1 



W2 + 



dH \ 



dH 



— (w5 + We) 



X, 



X„ 



W2 + 1^4 



q '^1 

Substituting Eqs. (Il2])-(l4l]) to Eqs. ([33]) and ([34]), we finally obtain 



0„ 



If U -V2\ If U-V2 

2 1--^ ^1 + 2 1 + -^) ^2 + ^3 



In the limit of g — )■ 1, Eqs. (1451) and ( 1461) reduce to 



1 - 
1 



If U-V2\ If U-V2 
W5 + we--\l+ ^ 1 wi--\l ^ ) W2-W4 



^1 



2 V D 



3ei 



U-V2 

D 



01) with 



(41) 



3^2 



— e 



In the limit of T = 0, the auto-correlation in the BGS and NES is given by 

U-V2 



D 



\ for {U - V2)/t = 0, 
1 for {U - V2)/t 00. 



(42) 
(43) 
(44) 

(45) 

(46) 

(47) 
(48) 

(49) 
(50) 



D. Concurrence 



The concurrence F has been proposed as a measure of entanglement for systems of two 
qubits [59I, l60| • It is defined with eigenvalues > ■ ■ > A4 for the positive Hermitean matrix 



R = VppVp by 



59 



60| 



r = max(Ai — A2 — A3 — A4, 0), 



(51) 



where p = {a^ a'^)p*{a^ ® cr^) and * denotes the complex conjugate. The entanglement of 
formation Ep 



801 1 is expressed in terms of F by 
Ep = 



5=± 



(52) 



For a pair of qubits of the Hubbard dimer, the concurrence is given by [IJ, |15| 



F = max ^1 ^ (^aj^a2<x^ | - ^ {nian2a) , 0^ 



(53) 



where the bracket (■) means the average over the density matrix. It is straightforward to 
generalize the method of Ref. [jj] to the NES, in which the g-dependent concurrence F^ is 
given by 



Tq = max 



max 



(Hf)j-(lf),4 



In deriving Eq. fl55|) . we employ the relations given by Eq. f l3T|) and by 

'dH_\ 

By using Eqs. ([5]) and (HTj) with 9 = t, v/e may calculate {dH/dt)q, 

/dH\ 8t 



(Wl - Uis)- 



Substituting Eqs. f HSj) and fl57j) to Eq. fl55j) . we obtain 



Fg = —max 



At 

— ]\ {wi-W2)\ - (W5 + We), 



In the limit of T = 0, F^ in both the NES and BGS is given by 

F =^ 

1 for (f/ - V2)/t = 0, 
for {U - V2)/t ^ oo. 
In the limit of q = 1.0, Eq. (l58ll reduces to 



^1 



-max 



(54) 
(55) 

(56) 

(57) 

(58) 

(59) 
(60) 

(61) 



which is in agreement with the result of Ref. ISj- With increasing the temperature, the 
concurrence is decreased and vanishes above the threshold temperature, as will be shown in 
Sec. 3B. 
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E. Conditional entropy 



The conditional entropy for subsystems A and B in the NES is expressed by 

S,iA,B)-S,iB) 



40 



45| 



with 



S,iA\B) 



Sg{A,B) = kB 



-il-q)S,iB)/kB-' 

Tr[p,(A,j?)]^-l 

1-q 
Tib [p,m' - 1 



S,{B) = kB 

1-q 

p,{B) = TiapM^B), 



(62) 



(63) 

(64) 
(65) 



where Tr^ stands for the partial trace over the state A and Pq{B) denotes the marginal 
density operator. In the limit of q — ?■ 1, 5*^(^15) reduces to the von Neumann conditional 
entropy, 5*1(^11?) = Si{A, B) — Si{B), whose properties have been discussed in Refs. 



61 



62|. 



In independent subsystems A and B w. 
dSD and dSD lead to Sg{A\B) = Sg{A) 



lere the relation: p{A, B) = p{A) ® p{B) holds, Eqs. 



40 



45|. 



Regarding subsystems A and B as sites 1 and 2 in the Hubbard dimer under consideration, 
we may obtain the marginal density operator given by 



p,{l) = Tr2p,(l,2) 

= ^ (^7i|0)i(0|i + g2\ t)i(t li + 93\ Did li + 9a\ tDiiU |i) 



with 



91 

92 
93 



94: = - {wi sin^ 9 + W2 cos^ 9 + w^) 
1 



[wi cos^ 9 + W2 sin^ 9 + 104 + 2^5) 
[wi cos^ 9 + W2 sin^ 9 + W4 + 2wq) 



(66) 

(67) 
(68) 
(69) 



where /3q(l, 2) = pg given by Eq. (^^. From Eqs. (I62|) -f l69|) . the conditional entropy is given 
by 



Rn 



Sg{l\2 



k 



B 



where 



(1-g) \dg 



d, = Tn [p,(l)]^ = Y^{9t + 9l + 9l + 9l) . 



(70) 

(71) 
(72) 



11 



When ei is the lowest eigenstate at T = 0.0, we obtain wi/Xg = 1.0 and Wi/Xg = 0.0 for 
z 7^ 1, which lead to 



p,(l,2) = 

f cos^e 
Pg(i) = f 



t)i(tii + u)iaii) + 



sin^^ 



i (|0)i(0|i + I t)i(t li + I Did li + I ti)i(ti li) for U/t ^ 0, 
HI t)i(t li + IDiUli) for f//t^oo. 



(73) 

(io)i(oii + in)i(nii), (74) 

(75) 



By using Eqs. (!70|) . ( !7T|) . (172|) and (!74|) . we obtain the conditional entropy given by 



Rn 



29-1 



(cos^'J 6' + sin 



2g, 



(4^-1 - 1) for U/t ^ 0, 
_kB_ (2?-i _ 1) for U/t oo. 



(76) 
(77) 



For q = 1.0 with t/U < 1.0 and \/i = 1/3 = 5 = 0, calculations using Eqs. (JUD and f lTSj) 
yield 



^ ^ ^ ^ 2t2 
Xi f/2 



(^1 _ Se-^'^*'/^) 



£ = f = ^1 - f^j (1 - 36-^'^/^ + 3e-^*^/^ 



(78) 
(79) 



from which the conditional entropy is given by 



with 



R, = 5i(l,2)-5i(l), 



5i(l,2) = 3ks[l + ^]e-'^''/^ 



(80) 



U 



SAl) = fcs ln2 + 



4^2 
LP 



1 - In ) - 6 e-^'^^V^ 



(81) 
(82) 



At T = 0.0, Eqs. 
quantum correlation [61 



yield -Ri = — A;Bln2 where the negative Ri expresses the 
62|. With raising the temperature, the conditional entropy is 
increased and changes its sign from negative to positive because of a contribution of the 
classical correlation, as will be shown is Sec. 3C. 
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III. MODEL CALCULATIONS 



A. Auto- and pair correlations 

We have performed numerical calculations, solving self-consistent equations for Eg and 
Xg given by Eqs. fl25|) - fl27j) . by using the Newton- Raphson method with Vi = V2 = B = 
otherwise noticed. Figure [1] shows the temperature dependences of various correlations 
given by Eqs. fl30|) - fl32|) and fl56|) for U/t = 5.0 with q = 0.6 and 1.0. With increasing the 
temperature, X]<T(^io-"'2o-)g is increased while both X]cr('^io-'^2-cr)g and X]cr('^L'^2o- + «L'^io-)g 
are decreased. In contrast, X]^=i(%t^j4)g is almost temperature independent. Temperature 
dependences of correlations for g = 0.6 are less significant than those for q = 1.0. These q 
and T dependences of correlations shown in Fig. [1] reflect on those of Og, Lg, Tg and Rg, as 
will be shown in the following. 

Figures [2] (a) and (b) show the temperature dependence of the auto-correlation (Og) for 
U/t = 0.0 and 5.0, respectively, with q = 0.6, 0.8, 1.0 and 1.2. The magnitude of spin 
correlation is given by (si ■ Ss), = (3/4)0,. At T = 0.0, Og is 0.5 and 0.89 for U/t = 0.0 
and 5.0, respectively, independently of q. Ogioi U/t = 0.0 is increased with increasing the 
temperature. For U/t = 5.0, Og is once increased with raising T, but it is decreased at 
higher temperatures after showing the maximum. 

Temperature dependences of pair correlation {—Lg) for U/t = 0.0 and 5.0 are shown in 
Figs. [2t^c) and (d), respectively, with q = 0.6, 0.8, 1.0 and 1.2. At T = 0.0, we obtain Lg = 
—0.5 and —0.89 for U/t = 0.0 and 5.0, respectively, independently of q: the negative sign of 
Lg stands for antiferromagnetic correlations for adopted parameters of Vi = V2 = B = 0.0. 
When the temperature is increased, magnitude of Lg is monotonously decreased. We note 
that —Lg for q < 1.0 is larger than that for q = 1.0 at ksT/t > 0.3, which expresses the 
intrigue effect of the nonextensivity on the pair correlation. 

B. Concurrence 

Figures [2t^e) and (f) show temperature dependences of the concurrence {Tg) for U /t = 0.0 
and 5.0, respectively, with q = 0.6, 0.8, 1.0 and 1.2. At T = 0.0, = 1.0 and 0.63 for 
U/t = 1.0 and 5.0, respectively, independently of q. With increasing the temperature, F, is 
more slowly decreased for smaller q. Tg vanishes above the threshold temperature Tp which 
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is implicitly determined by 

I ExpJ-A^(ei - Eg)] - Exp J- A, (62 - Eg)] \ 
-Exp^[-/3tHie, - Eg)] - ExpJ-A;,(e6 - Eg)] = 0, (83) 

with f3th = l/fcfiTr and Eg = Eq(Tr). Figures |2]^e) and (f) show that with decreasing q 
below unity, the threshold temperature Tr is increased. This is more clearly seen in Fig. [3] 
where the solid curve shows the q dependence of Tp for U/t = 5.0 (the dashed curve will be 
explained below). 



C. Conditional entropy 

Figures Wig) and (h) show temperature dependences of the conditional entropy (Rg) for 
U/t = 0.0 and 5.0, respectively, with various q values. Rg is negative at lows temperature 



which expresses the quantum entanglement 6l|], and it becomes positive at higher temper- 



atures. The threshold temperature Tr at which Rg changes its sign is implicitly expressed 
by 

CgiTR) = dgiTn). (84) 

The q dependence of Tr is plotted by the dashed curve in Fig. 3, which shows that Tr is 
increased with decreasing q below unity. Tr is correlated with Tp as shown in the inset of 
Fig. 4, although Tr does not agree with Tp. 

We note in Figs. [2]^a)-(h) that temperature dependences of O,, T^, and Rg become 
more significant with increasing g, which is consistent with the more significant temperature 
dependences in the specific heat (Cg) and susceptibility ixq) shown in Fig. El 



IV. DISCUSSION 



A. Relations among Lg, Tg and Rg 

It is interesting to investigate the relations among Lg, Vg and Rg. In Fig. Hl^a), Vg for 
U/t = 0.0 and 5.0 with g = 0.6 and 1.0 are plotted as a function of —Lg, which shows a linear 
relation: Tg ~ a{—Lg) —b{a,b>0). This linear relation between Tg and —Lg shown in Fig. 
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m^a) is realized in the parametric plot of —Lq{T) versus Tq{T) with fixed model parameters. 
However, it does not hold between —Lq and Tq when the model parameters are changed. 
This fact is easily realized when we compared Eqs. f lSOp with Eq. (160 p . which shows that 
with increasing U, \Lq\ is increased but Tq is decreased. In Fig. ID^b), Tq for U/t = 0.0 and 
5.0 with q = 0.6 and 1.0 is plotted as a function of —Rq, which shows the correlation between 
Tq and Rq. We note in Figs. HJ^a) and (b) that Lq, Tq and Rq are correlated although the 
precise relations among them are not clear. 



B. Effect of magnetic field and interatomic interactions 

We have so far assumed Vi = V2 = B = 0.0 for which the lowest eigenvalue of ei leads 
to the singlet ground state. If Vi, V2 and/or B are, however, introduced, the triplet state 
may be the ground state. The critical condition for the singlet-triplet transition is given by 
ei = 65, i.e., 

fiBB = ^(2Vi-V2-U+ ^{U - V2f + 16t2) . (85) 

Figure [5] shows temperature dependences of for q = 1.0 (dashed curves) and q = 0.6 
(solid curves) when B is changed with U/t = 5.0 and Vi = V2 = 0.0, for which Eq. (!85|) 
yields the critical field given by ^sBc/t = 0.351. The triplet state becomes the ground state 
for B > Be, where the pair correlation and marginal entropy are positive {Lq > 0, Rq > 0) 
and the concurrence vanishes (F^ = 0). 

Similarly, when we introduce Vi and/or V2 which satisfy Eq. fl85l) . the triplet state 
becomes the ground state even if S = 0. In the triplet state, we obtain Lg > 0, -R^ > and 
Tq = 0. The effect of interatomic interactions on energy, entropy and specific heat of the 



Hubbard dimer in the singlet state has been investigated within the NES 71 1. 



C. Physical meaning of the entropic index 



The entropic index q is conventionally regarded as a parameter which is determined by a 
fitting of the g-exponential distribution to experimental data except for some cases where q 
may be determined in a microscopic way 311] • We will briefiy discuss the physical meaning 
of the entropic index in a small system coupled to finite bath for which q is theoretically 



derived with the use of the MCA 



eikzl and ssA \mml 
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(1) MCA 

We consider a microcanonical ensemble consisting of a system and a hath with energies of 
Es and i^^, respectively {E = Es + Eb is conserved). Available states for the system with 
the energy between Es and Eg + ^Es are given by 



63 



65] 



(86) 



where fi^(i?) denotes the structure function expressing the number of states with energy E 
in ?7 (=S, B and S+B). We assume that the structure function is given by [63|, |65| 



Vt^{E) = Kmr,E"'^-\ 



(87) 



where K is a constant and the degree of freedom of variables in rj. Eq. (1871) is justified 
for d- dimensional A^-body ideal gases and harmonic oscillators, for which m = dN/2 and 
(iiV, respectively. For Es ^ Eb and nis <ti ttib, Eqs. (!86i) and (1871) lead to the PDF given 
by y, m] 



IP \ "T-B 



p{Es) « ( 1 - ^ 



a: Exp (-g/3E5), 



(88) 
(89) 



with 



<1 



1 + ^, 
rriB 

1 



(g-1)^' 

where Expg(a;) denotes the g-exponential function given by Eq. (12T1) . Eq. 
to the PDF obtained in the normal average. In the case of — > oo, Eq. 

p{Es) cx e-^""', 



(90) 
(91) 

(l89l) corresponds 
i8l) reduces to 

(92) 



with 



me 
E 



kBT 



(93) 



The specific heat of the bath is shown to be given by 6^ 



B 



dE_ 
16 



oc 



(94) 



Eqs. ([89]), ([90]), (192]) and dM]) imply that for finite bath, the PDF is given by the q- 
exponential function whereas for infinite bath (C„ = oo), the PDF is given by the exponential 
unction. A similar analysis has been made within the microcanonical approach in Refs. 



66 



67|. 



(2) SSA 

In the superstatistics 



pressed by a superposition over the spatial 



68|-|70|, it is assumed that properties of a given system may be ex- 
ly and /or temporarily fluctuating intensive pa- 
-|70|. Since the PDF of the equilibrium state i 



rameter (i.e., the inverse temperature) |68 
with the inverse temperature /3 (= l/ksT) is given by e~^'^'/Zi (/3), the PDF averaged over 



fluctuating /3- variable is assumed to be given by |68l470 1 
with 



(95) 



/(/3) 



n 



n/2 



(3 



n/2-lg-n/3/2/3o_ 



(96) 



r(n/2) V2/3o. 

Here T{x) is the gamma function and /(/3) denotes the x^-distribution with rank n which 
expresses the distribution of sum of squares of n random normal variables with zero mean and 
unit variance |69|. Average and variance of (3 are given by (/3)/3 = /3o and ((/3^)/3 — Pi)/ Pi = 



2/n. When adopting the type-A superstatistics in which the {3 dependence in 2'i(/3) is 
neglected [t^, we obtain (with i 



Pi oc 



1 H e, 

ns 



-ns/2 



(97) 



which is rewritten as 



Pi oc Exp (-g/3oei) 



(98) 



with 



1 - 



(99) 



Eq. (I98l) is in conformity with the normal- average PDF. 

It has been shown by a detailed microscopic calculation that the distribution of the inverse 
temperature of a system containing independent n particles coupled to a bath characterized 
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by a fixed inverse temperature of /3, is given by 81 1 



r(n/2) V 2 

Eq. f llOOp expresses the inverse-gamma distribution, and its profile is similar to that of the 
gamma-distribution given by Eq. ( !96|) for large n [sij. Unfortunately we cannot obtain an 
analytic expression for the PDF averaged over /t(/3) by using Eq. with — )■ /t(/9). 



Nevertheless the calculation of Ref. 



V. CONCLUSION 



8l| justifies the concept of the superstatistics. 



We have calculated various quantities of quantum entanglement such as auto- (Og) and 
pair correlations (Lg), concurrence (Tg) and conditional entropy (Rq) of the half- filled Hub- 
bard dimer as functions of the entropic index and the temperature within the framework of 



the NES [28[. It has been shown that the properties of Og, Lg, Tg and Rq are considerably 
modified by the nonextensivity. In particular, for q < 1.0, the thermal entanglement my be 
survive at higher temperatures than that for q = 1.0, because the threshold temperature 
where Tg vanishes (Tr) or Rq changes its sign (Tr) is more raised for a larger (1 — q) (Fig. 
3). The three measures of Lg, Tg and Rg for thermal entanglement are correlated each other 
although the precise relations among them are not clear. 

The NES has a wider applicability than the BGS, which corresponds to the q = 1.0 limit 
of the NES. We note that the PDF in the MCA given by Eq. (jSH]) is equivalent to that in 

There is, however, distinct differences in their expressions of q 



[Eqs. (EO]) and (ES])] 



s. ( 


981) 


82 


-84 



1 + ^ > 1.0 in the MCA, 

(101) 

1 - ^ < 1.0 in the SSA. 

ns — 

The entropic index in the MCA is expressed in terms of the bulk parameter (m^) while that 
in the SSA is given in terms of the system one (ns). Furthermore the conceivable value of q 
in the MCA is different from that in the SSA. In this respect we have not obtained a unified 
physical interpretation of the entropic index at the moment. Nevertheless Eq. fllOll) shows 
that the entropic index q may be related with the size of the system and/or bath and that 
the nonextensivity may be realized in such a small-scale system. It might be interesting to 



18 



perform experiments by changing the size of the system and/or bath, in order to examine 
the possibihty that the nonextensivity reflects on the thermal entanglement of two-qubit 
Hubbard dimer. Such experimental studies might clarify the role of the nonextensivity in 
small systems and provide valuable insight on the validity of the MCA and SSA. 
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Appendix: A. Thermodynamical quantities of the Hubbard dimers with the original 
MEM 

1. Energy and entropy 

When the energy Eg and partition function Xg are obtained by solving Eqs. f l2^ - fl27]) 
for given q and /3, the Tsallis entropy given by Eq. ([T]) may be calculated by 

S, = kB( ^1 , (A.l) 



1-q 



with 



2. Specific heat 



The specific heat is given by [55|, |56| 



dl3\ fdEg 



where dEq/d(3 may be determined by simultaneous equations given by 

dE„ dEg dXg ^ 

dXg dEg dXg ^ 

—S. = a2i—^ + a22-T-^ + b2, (A. 5) 

dp dp ' ^ ' 
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with 

an = y^^j^"'' Ci, ai2 = ^21 = PXq, 022 = 0, 

h = -^J^^r'^e.le.-E,), 62 = 0. (A.6) 



Eqs. (lAa-lira are derived from Eqs. (^-^. Solving Eqs. (lA^-lfAlHl) for dEg/df3, we 



obtain 



KbI J i- - an - cn2a2i 



In the hmit of g = 1.0, Cg becomes 

n. = I 

3. Susceptibility 



Ci = ("7^2 (A.8) 



The paramagnetic spin susceptibihty Xg is given by |55l. l56| 



where = d'^EjdB^\B=o and xf^ = S^Xg/SE^Is^o- From Eqs. i^E)-<^, we obtain 
the simultaneous equations for Eg^^ = d'^Eq/dB'^\B=o and Xg"^^ = d'^Xg/dB^\B=o, 

= aiiEf + ai2Xf ) + A, (A.IO) 
= a2ii?f + a22Xf + /2, (A.ll) 



with 



f2 = -Yl (A.12) 

From Eqs. flX9l) - f[02l) . we obtain 

021 ~^ 

which does not include j\. In the limit of g = 1.0, the spin susceptibility is given by 
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4. Model calculations 



Fig. [6] shows temperature dependences of the calculated entropy, specific heat and sus- 
ceptibility. Temperature dependences of the entropy 5*^ for U /t = 0.0 and 5.0 are plotted in 
Figs. 6(a) and (d), respectively, with q = 0.6, 0.8, 1.0 and 1.2. For a larger U value, Sg is 
more quickly increased at low temperatures. This is because the energy difference between 
the ground state (ei) and the first-excited state (64) becomes smaller when the strength of 
U is more increased. With increasing q, the saturation value of Sq at higher temperatures 
becomes smaller. Temperature dependences of specific heat Cq for U/t = 0.0 and 5.0 are 
plotted in Figs. 6(b) and (e), respectively, for various q values. Figs. 6(c) and (f) show 
temperature dependences of the susceptibility Xg ^oi U/t = 0.0 and 5.0, respectively, with 
q = 0.6, 0.8, 1.0 and 1.2. We note that temperature dependences of Cq and Xg at low tem- 
peratures for q = 1.2 are more significant than those of g = 1.0 whereas those for g = 0.8 is 
less significant than those of g = 1.0: temperature dependences of Cq and Xg become more 
significant with increasing q. In contrast, Sg is increased with increasing q. These behaviors 
are understood as follows. The q- and T-dependent thermodynamical quantity Qq{T) may 
be expand at g = 1.0, 

^ Q,{T) + {q-l)Q^^\T) + -, (A.16) 

where Qf^\T) = d^Qq{T) / dq^\q=i. Actual analytical evaluation of Q^^\t) is tedious be- 
cause it involves self-consistent calculations as discussed in preceding subsections. Our model 
calculations show that Q'^i\t) > for Cq and Xg whereas Q^i\t) < for Sq at low temper- 
atures. The characteristic temperature dependences in thermodynamical quantities depend 
on the entropic index g. When comparing these results with the counterparts obtained in 



our previous study [55] using the normalized MEM with g-average 30| , we realize that both 



results approximately have the g ^ 1/g symmetry: for ex ample, results of g = 0.6 in Fig. [6] 



are similar to those of g = 1.5 in Figs. 2, 3 and 4 of Ref. 55|, |79 |. 
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Appendix: B. Heisenberg dimers with the original MEM 



It is well known that the Hubbard dimer with U/t^ 1 (with Vi = V2 = 0) is equivalent 
to the Heisenberg dimer with the superexchan ge in teraction of J^e ~ —t^/U. It is worthwhile 



to apply the NES with the original MEM 



28 



73| to a Heisenberg dimer given by (s = 1/2) 



H = -Jsi-S2-gfiBB{si^ + S2z), (Bl) 

where the exchange interaction J is positive (negative) for ferromagnetic (antiferromagnetic) 
coupling, g (=2) denotes the g-factor, fiB the Bohr magneton, and B an applied magnetic 
field. Four eigenvalues of H are given by 

ei = —-q^bB, ~ + g^BB for i = 1, 2, 3 and 4, (B2) 

and corresponding eigenvectors are given by 

i)i|t)2), (B3) 

i)i|t)2), (B4) 

(B5) 
(B6) 



$1) = 


^(lt).U)2 


$2) = 


;^(lt).|i)2 


$3) = 


1 t)l|t)2, 


$4) = 


li)l|i)2. 



It)2- 




i)2) 


1 t)2 + 




\)2) 



(B7) 
(B8) 



When 5 = 0, |$3) and |$4) may be alternatively expressed by 

i*=> ^ ^(1 
I*.) ^ ^(1 

Eqs. dHU), dHl]), dHl]) and ([BS]) form the Bell basis. 
In the BGS the partition function is given by 

Zbg = Z^ = e-^ + e^ [1 + 2cosh (2/3/xbB)] , (B9) 

from which various thermodynamical quantities are easily calculated. The susceptibility is, 
for example, given by 
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The calculation for the Heisenberg dimer within the NES goes parallel to that for the 
Hubbard dimer with the use of four eigenvalues given by Eq. flB2l) . The average energy and 
partition function are expressed by 



(Bll) 
(B12) 



with 



ExpJ-/3(e,-E,)]. 



(B13) 



The pair-correlation function (Lg) is defined by 



-E 



wl - 



dJ 



which yields (with 5 = 0.0 hereafter) 



^{(Exp 



Exp„ 



(B14) 
(B15) 
(B16) 



For g = 1.0 with 5 = 0, Eq. (]B17p reduces to 



/ £J 313 J \ 

( e 4 — e 4 1 . 



1 / M 3/3 J' 

Li = — ( e 4 - e 4 



}. (B17) 



(B18) 



The concurrence for J < is given by 



Tq = —max 



2 wi - ^ Wi\,0 



(B19) 



For g = 1.0 and B = 0, Eq. ( 1B19I) becomes 



r„ = — max [| 



3/3 J ;3J 

e 4 - 3e 4 



(B20) 



which is in agreement with the result of Ref. 
By using the marginal density matrix given by 



P.(l) = 2(lt)i(t|i + li)i(i|i) 



(B21) 
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we obtain the conditional entropy [40| 



Rr, 



■(29-1X1-" - 1) 
. (1-9) 



(B22) 
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FIG. 1: (Color online) (a) Temperature dependences of various correlations for U/t = 5.0 with 
q = 0.6 (solid curves) and 1.0 (dashed curves). 

FIG. 2: (Color online) Temperature dependences of auto-correlation (Og) for (a) U/t = 0.0 and 
(b) U/t = 5.0, pair correlation {—Lq) for (c) U/t = 0.0, (d) U/t = 5.0, concurrence (Tq) for (e) 
U/t = 0.0 and (f) U/t = 5.0, and the conditional entropy {Rq) for (g) U/t = 0.0 and (h) U/t = 5.0 
with q = 0.6 (solid curves), 0.8 (dotted curves), 1.0 (dashed curves) and 1.2 (chain curves). 

FIG. 4: (Color online) (a) Fg as a function of —Lq and (b) Fq as a function of —Rq with q = 
0.6 (solid curves) and 1.0 (dashed curves), the result oi U/t = 5.0 in (a) for q = 0.6 being 
indistinguishable from that for q = 1.0. 



FIG. 5: (Color online) Temperature dependences of the pair correlation Lq for various B values 
for U/t = 5.0 with q = 0.6 (solid curves) and 1.0 (dashed curves). 



FIG. 6: (Color online) The temperature dependences of (a) Sq, (b) Cq and (c) Xq fo^ U/t = 0, and 
(d) Sq, (e) Cq and (f) Xq fo^ U/t = 5 with q = 0.6 (solid curves), 0.8 (dotted curves), 1.0 (dashed 
curves) and 1.2 (chain curves). 



FIG. 3: (Color online) q dependences of threshold temperatures of Tr (the solid curve) and Tr 
(the dashed curve) for U/t = 5.0, the inset showing Tr against Tp. 
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